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The  Chemistry  of  Detonations.  1^.  A  Simple  Method  for  Calculating 
Detonation  Properties  of  C-H-N-0  Explosives 


Mortimer  J.  Kamlet  and  S.  J.  Jacobs 
U.  S.  Naval  Ordnance  Laboratory,  White  Oak,  Silver  Spring,  Maryland 

ABSTRACT.-  Detonation  pressures  of  C-H-N-0  explosives  at  initial  densities 

above  1.0  g/ce  may  be  calculated  by  means  of  the  simple  empirical  equation, 

2  1  X 

P  =  Kp  cp  ,  K  =  15.58,  9  =  NMZ  Q2  |  detonation  velocities  by  the  equation,  D  = 
A  0 

A©2  (1  +  Bp  ),  A  =  1.01,  B  =  1,30.  N  is  the  number  of  moles  of  gaseous  detona- 
0 

tion  products  per  gram  of  explosive}  M  is  the  average  molecular  weight  of  these 

gases}  Q  is  the  chemical  energy  of  the  detonation  reaction  (AH0  per  gram)}  and 

p  is  the  initial  density.  Values  of  N,  M  and  Q  may  be  estimated  from  the 
0 

[Ha0-C0aJ  arbitrary  decomposition  assunmtion,  so  that  the  calculations  require 
no  other  input  information  than  the  explosive's  elemental  c  cm  position,  heat  of 
formation  and  loading  density.  Detonation  pressures  derived  in  this  manner 
correspond  quits  closely  to  values  predicted  by  a  computer  code  known  as  RDBT, 
which  employe  the  most  recent  parameters  and  co volume  factors  with  the 
Kiatiakowaky-Wilscn  equation  of  state. 
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THE  CHEMISTRY  OF  DETONATIONS.  I.  A  SIMPLE  METHOD  FOR  CALCULATING 
DETONATION  PROPERTIES  OF  C-H-N-0  EXPLOSIVES. 

This  report  is  the  first  of  a  series  which  describes  simplified 
methods  of  predicting  detonation  parameters  and,  eventually,  certain 
types  of  damage  effects  of  C-H-N-0  high  explosives  using  as  a  priori 
information  only  a  knowledge  of  their  chemical  structure.  It  is 
hoped  that  the  relationships  described  herein  will  be  useful  to  the 
synthesis  chemist  in  designing  new  more -efficient  explosives.  The 
work  was  carried  out  under  the  Foundational  Research  Program  of  this 
Laboratory. 
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I.  INl'RODUCTION 

To  a  chemist  concerned  with  the  synthesis  of  new  high  explosive  compounds 
the  ability  to  compute  detonation  properties  (detonation  pressure,  energy*  and 
velocity  as  well  as  product  composition)  frcm  a  given  molecular  structure  ani 
the  known  or  estimated  crystal  density  is  a  problem  of  the  utmost  importance* 

The  calculated  properties  could  be  meaningful  in  the  decision  as  to  whether  it 
is  worth  the  effort  to  attempt  a  new  and  complex  synthesis.  One  reason  behind 
the  recent  development  of  detonation  properties  programs  for  use  on  high  speed 
computers  has  been  to  supply  this  desired  information.  One  such  program,  the 
RUBY  code1,  ha 3  recently  been  made  available  to  a  number  of  laboratories,  the 
authors’  included. 

In  an  effort  to  understand  the  formidable  appearing  output  of  many  compu¬ 
tations  for  a  wide  variety  of  C-H-H-0  explosives  at  various  initial  loading 
densities,  we  have  investigated  interrelationships  between  such  properties  as 
pressure,  velocity,  density,  heat  of  reaction,  etc.  These  studies  have  led  to 
a  number  of  interesting  observations,  important  among  which  were  the  "facts’’  that 
much  simpler  semi-empirical  formulas  could  be  written  for  desk  calculation  of 
detonation  velocities  and  detonation  pressures,  with  about  the  same  reliance  an 
their  answers  as  one  could  attach  to  the  more  complex  computer  output.  These 
equations  require  as  input  information  only  the  explosive's  composition  and 
loading  density  and  an  estimate  of  its  heat  of  formation  and,  in  their  oonroarative 
simplicity,  seem  to  thrcv  light  on  the  relative  importance  of  the  quantities 
which  determine  the  detonation  pressure  in  particular,  and  other  properties  as 
well. 

It  ie  hoped  that  the  present  findings  may  give  those  untrained  in  the 
details  of  the  Uieraodynamic-hydrodynamic  calculations  a  better  "feel”  far  the 
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results  of  the  computer  output.  A  further  hope  is  that  these  papers  will  also 
serve  as  a  reminder  to  users  of  11  black  box"  computation  schemes  to  the  effect 
that  "the  output  is  no  bettor  than  the  input  and  the  responsibility  for  the  input 
still  rests  largely  with  the  user,  not  the  box  maker". 

II.  DETONATION  CALCULATION 

In  the  last  25  years,  calculations  of  {ho  detonation  properties  of  con- 
densed  explosives  from  thoir  chemical  compositions  and  densities  have  been 

3 

approached  in  various  ways.  All  have  used  the  necessary  conservation  conditions 
for  steady  flow  with  the  detonation  discontinuity  satisfying  the  Chapman-Jouguet 
hypothesis  (minimum  detonation  velocity  compatible  with  the  conservation  condi¬ 
tions  or  sonic  flow  behind  the  discontinuity  in  a  reference  frame  where  the 
discontinuity  is  at  rest).  In  order  to  describe  the  product  state  and  the  thermo¬ 
dynamic  variables  whic.„  fix  its  composition,  an  equation  of  state  applicable  to 
a  vary  dense  state  is  required.  To  apply  this  equation  to  a  mixture  of  gaseous 
and  solid  products,  «  mixing  rule  is  also  needed  and  the  temperature  must  be 
explicitly  defined.  Consequently,  the  choice  of  equation  of  state  to  be  used  for 
the  "gaseous"  products  must  be  somewhat  more  general  than  the  equations  used  to 
describe  properties  on  an  isentrope  from  the  detonation  state.  Of  the  several 
equation  of  state  approaches  used,  only  three  will  be  mentioned  here. 

First,  wo  consider  the  virial  expansion  in  density  originally  due  to 
Boltzmann*’5  derived  from  the  kinetic  theory  of  gases  for  hard  sphere  molecules. 
This  equation  was  modified  by  Hirschfeldor  and  Rosevoare5  and  covolume  terms  for 
product  species  were  adjusted  to  high  temperature4*5  by  setting  them  equal  to 
the  high  temperature  second  virial  coefficients.  The  equation  is: 

FV  /RT  =  1  +  x  +  .625  x2  +  .2869  x3  +  .1928  x4  (la) 
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x  =  b/V  >  b  =  Ix1bi  (lb) 

5.3  the  molar  covolumo  of  the  ith  specie, 

x^  is  the  mols  fraction  of  the  ith  specie, 

V  is  the  molar  volume  of  the  gas  mixture,  and 
E 

P,  T  and  R  are  the  pressure,  temperature,  and  gas  constant*. 

The  b^  are  derived  from  the  collision  radii  of  the  molecular  3pocie3  at 
high  temperature  and,  as  in  the  kinotic  theory  of  gases  at  moderate  pressure,  are 
equal  to  four  times  the  molecular  volume  multiplied  by  Avogadro's  number.  Deapite 
the  use  of  diminished  covolumes  in  the  equation  and  despite  the  aonarent  theore¬ 
tical  basis  of  the  model,  the  equation  is  oversimplified  and  the  results  on 
detonation  calculations  quite  clearly  show  it  to  be  inaccurate. 

The  virial  expansion  in  the  pressure  was  used  by  Jones  and  Miller6. 

Their  equation  is; 

W  =  RT  +  bP  +  cP2  +  dP3  .  (2) 

g 

The  constants  b,  c,  and  d  for  this  equation  wore  fitted  by  the  authors  to  give 

the  correct  detonation  velocity  (D)  vs  initial  density  (p  )  behavior  for  the 

0 

explosive  described  (TNT).  Although  this  relationship  is  empirical  and  the 
constants  are  fixed  by  detonation  data,  the  equation  exhibits  solid-like  proper¬ 
ties,  e.  g.,  an  internal  oner©.'  term  due  to  molecular  repulsion  at  high  pressure, 
and  to  this  extent  appears  to  describe  the  state  of  a  high  density  gas  better 
than  Equation  (l).  That  this  equation  has  not  been  widely  used  by  others  may 
have  as  its  reason  that  the  constants  may  not  bo  germ  rally  applicable  to  all 
compositions.  If  the  constants  require  known  detonation  velocity  data  for  their 
determination,  the  equation  would  not  be  suitable  for  ji  rriorj.  calculations  from 
composition  and  density. 

A  third  approach,  instigated  by  Kistiakcwsky,  Wilson  and  Halverson7,  nay 
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be  said  to  hava  its  roots  in  Eq,  (l).  These  author's  modified  an  equation  due 
to  Becker, 

P  =  RT(l  +XeX)/V  +  f  (V)  1  X  =  t/V,  (3) 

by  dropping  the  f(V)  dependence,  adding  an  adjustable  constant,  0,  and  making 
b  a  function  of  temperature.  This  "variable  covolume"  equation,  as  further 
modified  by  Fickstt  and  Cowan0 ,  became: 

FVg/RT  =  1  +Xe(JX  ,  (4a) 

X  =  Klx^/Vgd  +  8)a  .  (4b) 

Eq.  (4)  is  a  variable  cc volume  departure  from  the  hard-sphere-molecule  Eq.  (1), 
for  if  p  =  0.625,  K  '  1  and  a  =  0  the  K-W  equation  would  be  identical  with  the 
Boltzmann  equation  to  the  third  viriai  term  and  the  k^s  would  be  Just  the  b^a 
of  tiie  hard-sphere-molecule  model.  If  p  were  0.625,  with  a  about  0.25  to  0.5, 
one  might  consider  the  K-W  equation  to  be  a  "soft>-3phercn  equation  of  state.  In 
applying  Eq.  (4)  to  calculation  of  detonation  velocities  it  was  quickly  found, 
however,  that  p  could  not  be  as  large  as  0.625  and,  in  the  earlier  paper's  on 
this  problem7*9*10,  the  values  adopted  were  k  =  1.0,  p  =  0.3,  a  =  0.25  aid  8=0. 
With  these  parameters  arvd  the  D-p  data  for  a  number  of  explosives,  k.  's 

0  1 

were  determined  for  the  principal  molecular  species  expected  as  detonation 
products  from  C-H-N-0  explosives.  The  values  obtained  came  fairly  cloee  to 
agreeing  with  the  b^'s  of  Eq.  (l)  if  one  defined  b^  as, 

\  =  k  ^/(T  +  6)a,  (5) 

and  assumed  T  =  4000°  as  typical  of  detonation  temperatures  found.  Computations 
were  made  by  estimating  fixed  detonation  product  compositions  as  well  as  on  the 
basis  of  equilibrium  calculations.  The  equilibrium  calculations  of  Brinkley  and 
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Wilson9  tended,  at  that  time,  to  favor  an  [Ha0-C0-C0a]  "arbitrary"  method  of 
estimating  detonation  product  compositions.  Consequently,  Snay  and  Christian10 
tried  both  [Ha0-CQ-C0a]  and  [C0-Ha0-C0a]  arbitrary  decomposition  schemes  to 
test  the  affects  of  changing  composition  on  predicted  detonation  properties  at-i 
covolume  factors.  By  least  squares  they  determined  a  beet  set  of  covolume  fac¬ 
tors  for  the  above  parameters.  The  results  were  not  very  much  influenced  by 
the  decomposition  assumption,  but  computed  detonation  pressures  were  lower  than 
values  found  experimentally. 

Tne  next  step  toward  better  fitting  of  the  K-W  equation  to  detonation 
data  was  mace  by  Cowan  and  Pickett®,  who  established  a  substantially  different 
set  of  parameters  and  covolume  factors.  More  recent  adjustments  by  Mader11 
have  led  to  the  parameter  sets  used  most  frequently  tcday  in  the  RUBI  code. 
Mader's  parameters  were  designed  to  give  the  best  match  with  five  experimental 
measurements  considered  to  be  highly  accurate;  the  detonation  pressure  of  RDX  at 
1.3  g/cc,  the  detaration  velocities  of  RDX  at  1.0  and  1.8  g/cc,  and  the  detona¬ 
tion  velocities  of  TNT  at  1.0  and  1.64  g/cc.  Fundamental  difficulties  in  finding 
a  single  set  of  parameters  to  accomodate  these  five  measurements  led  Mader  to 
suggest  dual  sets  of  (3  and  k:  an  "RDX  parameter  set"  to  be  used  with  com  pounds 
producing  lesser  amounts  of  solid  carbon  in  the  detonation,  a  "TNT  parameter 
set"  with  explosives  producing  greater  amounts  of  solid  carbon  (Table  I). 


Table  I.  Parameters  in  hi.  (4) 


Source 

ref . 

a 

P- 

Y 

e 

Kistiakc-^^y-Wilsan  7*1 

Brlnkley-Wil 9 an 

9r 

0.25 

0.3 

1.0 

0 

Chris tian-Snoy 

10J 

Cowaa-Fickett 

8 

0.5 

0.09 

11.85 

400 

Mader,  RDX 

11 

0.5 

0.16 

10.91 

400 

Mader,  TNT 

11 

0.5 

0.096 

12.69 

400 

5 


NOLTR  67-66 


It  ohoula  be  said  that  only  since  Cowan  and  Fickett’s  report  has  good 
experimental  detonation  pro car re  data  been  available  for  use  in  arriving  at 
best-fit  parameters.  Al«o,  the  covolume  factors  (k, )  have  been  deliberately 
normalized  but  not  set  equal  to  molecular  excluded  volumes.  In  this  normaliza¬ 
tion,  rather  coger.t  arguments  were  made  for  reducing  the  orientation  effect  of 
the  polar  molecules  ( Ha0  and  NH,  in  particular) ,  and  thereby  increasing  the 
magnitudes  of  the  corresponding  kJs  relative  to  those  for  non-polar  molecules1*. 
The  result  is  an  timpird  equation  which  at  this  time  is  the  best  available 
for  general  detonation  state  calculations. 

Three  significant  consequences  in  regard  to  detonation  calculations  on 
C-H-N-0  compositions  derive  from  the  present  K-W  parameters  and  covolume  factors: 
the  predicted  detonation  temperature  iu  quite  lew;  the  predominant  carbon-oxygen 
product  is  C0a  rather  than  CO  over  a  vide  range  of  compositions  at  the  higher 
loading  densities;  the  detonation  pressure  ana  velocity  are  reasonably  close  to 
experiment  vher  experimental  data  are  available.  Whether  the  result  on  C0a  is 
correct  or  nw,  it  has  served  as  one  point  of  departure  for  the  discussion  which 
follow 8 • 

III.  THE  CQVOLUKE  FACTOR  -  MOLECULAR  WEIGHT  RELATIONSHIP 

The  covolumo  factors  most  recently  used  by  liader1 1  for  Ha0,  C0a  and  Na 
(and  thus,  by  inference,  for  isoelectronic  CO)  wore  further  adjusted  from  Cowan 
and  Fickett's  values8  so  as  best  to  reproduce  experimental  Hugoniots13”15,  i.  e. 
Ha0,  360  250;  C0a,  670  ->  600;  N„  330  -*■  380;  CO,  390  ->  390.  The  ki's  for  the 

"minor"  detonation  specie;:,  CH*,  Ha,  NO  and  0a,  in  current  K-W  ue lunation  ccmpu- 
tations  remain  the  "theoretical"  values  deriving  from  calculated  molecular  dimen¬ 
sions.  In  the  light  op  this  fact  and  because  Ha0,  CGa,  Na  and  CO  are  usually 
eonaidorad  to  comprise  %+%  of  the  detonation  gases  from  organic  C-Il-N-0  exploaivea 
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it  is  noteworthy  that  Mader's  k^'s  for  ths  "major"  detonation  specie®  are  closely 
proportional  to  thoir  molecular  weights  (Table  II). 


Table  II.  The  Covolume  Factor-Molecular  Weight 
Relationship  for  the  Major  Detonation  Species 


Specio 

-  _ki 

_____  Vi 

Ha0 

250 

13.89 

380 

'3.57 

CO 

390 

13.93 

coa 

600 

13.64 

average 

(13.76  +  0, 

This  observation  is  of  substantial  interest  in  view  of  Christian  and  Snay’s 
report10  that  empirical  covolume  factors  for  the  total  gas  mixtures  of  twenty- 
organic  explosives,  chosen  to  give  best  average  agreement  between  calculated  and 
measured  detonation  velocities,  were  also  roughly  proportional  (+  ca  4$)  to  the 
average  molecular  weights  of  the  assumed  gaseous  detonation  products.  Taken  in 
combination,  these  findings  have  suggested  that,  for  most  organic  C-H-N-0  explo¬ 
sives,  Lx^k^  in  Kistiak  iwsky-Wilson  equation  may  be  replaced  by  the  product 
H*M,  where  H  ib  the  constant,  13.7b,  and  M  is  the  average  ga s  molecular  wei^it. 

Since  (m/V  )  =  p  ,  p  being  the  density  of  the  gaseous  products  in  the  detonation 
S  S  £> 

state,  this  allows  transformation  of  Eq.  (4)  to: 


P  =  (RTp  /M)(1  +X«PX),  X  =  wHp  /  (T  +  6)°  . 

6  o 


(6) 


We  now  introduce  two  terms  which  will  become  very  important  in  subsequent 
discussions:  N,  the  number  of  moles  of  gaseous  detonation  products  per  gram  of 
explosive,  and  G,  the  weight  fraction  of  explosive  going  over  to  gaseous  products. 
From  the  definitions,  N*H  -  g,  so  that  after  combining  constants,  tcH  =  A  and 
{3kH  -  B,  the  equation  may  now  take  the  gonoral  form: 
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P 


N  R  T  p  A  p 

_ JL-  r*  +  - _ il, 

G  L 1  (T  +  6)a 


6Bp8/(t  +  a)0-, 


(7) 


Tho  various  factors  influencing  P  in  the  K-W  equation  may  be  a33eaced 
more  readily  when  it  takes  the  form  of  Eq.  (7)  in  which  the  k^  terns  are  elimi¬ 
nated  than  had  been  the  case  with  Eq<  (4).  As  an  example,  the  dependence  of 
P  on  T  resulting  from  the  use  of  Kador's  "RDX  parameters”  is  evaluated  at  various 
loading  densities  in  Appendix  I  for  a  special  case  of  Eq.  (7). 

IV.  THE  RUBY  COMPUTER  CODE 

With  the  advent  of  the  new  high-speed  digital  computers,  several  more 
sophisticated  methods  of  predicting  detonation  properties  of  high  explosives  fran 
the  Kistiakowsky-Wilson  equation  of  state  have  achieved  wide  acceptance.  Two 
related  programs,  the  STRETCH  BKW  code  for  use  on  the  IBM-7030  computer11  and 
the  RUBY  code1  for  use  on  the  IBM-7090  provide  the  comparison  information  against 
which  predictions  from  the  equations  offered  in  the  present  paper  will  be  judged. 
For  the  purposes  of  present  discussions,  RUBY  results  as  determined  at  the  Navrl 
Ordnance  Laboratory16  and  STRETCH  BKW  results  differ  only  in  minor  regards17. 
Unless  otherwise  specified,  therefore,  the  term  RUBY  shall  hereafter  encompass 
the  results  of  both  systems  and  RUBY  computations  shall  be  considered  as  based 
on  Mader's  most  recent  parameters  and  covolume  factors  in  Eq.  (4) ,  with  the  heat 
of  formation  of  solid  carbon  taken  as  zero. 

The  RUBY  code  finds  "exact"  equilibrium  compositions  of  the  detonation 
products  ty  sophisticated  multi-iterative  processes  involving  minimization  of 
tho  Gibbs  free  energy  of  the  total  system.  Complex  input  equations,  derived 
through  the  equation  of  state,  relate  free  energies  of  formation  of  the  potential 
detonation  products  to  temperatures,  pressures,  ccvolumes  and  the  compressibility 
of  solid  carbon.  The  current  codes  are  equipped  to  consider  as  many  as  sixteen 
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gasaoua  products  and  two  solid  phaso3;  improvod  codes  are  being  designed  to 
accanodate  sixty  gaseous  species  and  fivo  solid  phacos,a. 

Several  pages  of  print-out  are  associated  with  a  single  RUBY  computation. 
Reported  together  with  D,  P,  T  and  V  at  a  given  p  (loading  density)  are  the 
following  quantities:  AE0,  the  chemical  energy  of  the  detonation  reaction  (*»  Q)| 

Ej  -  E0,  the  change  in  internal  energy  across  the  detonation  front;  pj,  the 
Chapnan-Jouguet  density;  N,  the  total  number  of  moles  of  gaseous  detonation 
products  per  gram  of  explosive ;  N4 ,  the  number  of  moles  of  the  individual  gas 
species;  N  ,  the  number  of  gram-atoms  of  solid  carbon  or  other  solids;  V  ,  the 

S  3 

solid  volume  per  gram-atom  carbon;  and  y,  the  "gamma  law"  constant,  y  = 

(-  b  In  P/  6  In  V)  .  Although  only  P,  D  and  possibly  T  are  subject  to  relatively 
unambiguous  experimental  verification  at  the  current  "state-of-the-art",  RUBY1 a 
predictions  of  other  of  tho  above  quantities  provide  a  framework  upon  which  ouch 
of  the  subsequent  discussion  is  based. 

V.  SIMPLE  EQUATIONS  FOR  ESTIMATING  DETONATION  PROPERTIES 
We  have  found  that  estimates  of  detonation  pressure  and  velocity,  corre¬ 
sponding  surprisingly  closely  to  RUBY  predictions,  are  possible  for  C-H-N-0 
explosives  by  means  of  relatively  simple  empirical  equations.  These  equations 
imply  that  the  "mechanical"  properties  of  the  detonation  depend  only  on  the  number 
of  moles  of  detonation  gases  per  unit  weight  explosive,  the  average  molecular 
weight  of  these  gases,  the  chemical  energy  of  the  detonation  reaction  (Q  =  AH0)» 
and  the  loading  density,  with  the  dependence  in  each  case  being  relatively  simple. 
The  equations  take  the  forms, 

F  =  Kp29  ,  X=  15.58,  9  =  NM^Q^  ,  (8) 

0 

and 

=  +  A -1.01,  B  =  1 .30 

0 
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whore  P  is  expressed  in  Kbar,  D  in  mr/usec,  N  in  moles  g&a/ g  explosive,  M  in 

g  gas/ mole  gas,  Q  in  c&l/g  and  p  in  g/cc. 

0 

As  a  first  test  of  these  equations,  values  of  N,  M  and  Q  are  taken  from 
RUBY  print-outs  as  reported  by  Hurwitz16  to  compare  (EA«  E)  an^ 

^calc  (®9*  ^  Pjjuby  an^  ^RUBY  ^or  S0Vara^  representative  explosives  at 
typical  leading  densities.  The  results  are  given  in  Table  III. 

Table  III  near  here. 

It  is  seen  that  substituting  NRupy»  arid  into  E<13,  (E)  and  (9) 

leads  to  detonation  pressures  and  velocities  which  differ  only  nominally  from 
values  predicted  by  the  computer.  The  problem  of  a  suitable  method  for  simple 
hand-calculation  of  these  detonation  properties  resolves  itself  to  finding  a 
reasonable  scheme  for  estimating  N,  M  and  Q  without  the  assistance  of  the  computer 

VI.  ESTIMATION  OF  N,  M  AND  Qj  THE  [Ha0-C0*]  ARBITRARY 
Product  compositions  at  the  Chapman-Jouguet  state  and  in  the  subsequent 
expansion  of  the  detonation  gases  depend  most  strongly  on  the  two  important 
equilibria, 


2  CO  r— 

•  CO^  C  f 

AH0  =-41.2  Kcal, 

(10) 

Ha  +  CO  — ; 

— »  Ha0  +  C  , 

AH0  =01.4  Kcal. 

(ID 

Fran  Le  Cha telier's  principle,  higher  pressures  (higher  densities)  should  shift 
these  equilibria  to  the  right,  higher  temperatures  to  the  left.  Since  RUBY  treats 
carbon  as  a  condensed  phase  (and  thus  of  unit  activity; ,  it3  amount  does  not 
materially  affect  the  equilibria  and,  30  long  as  at  least  same  solid  carbon 
appears,  the  ratios  (C02/C0)  and  (HaO/Ha)  are  rough  measures  of  the  equilibrium 
positions.  RUBY's  predictions  of  these  ratios  for  acme  typical  explosives  are 
given  in  Tahle  IV. 
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Table  IV.  Detonation  Product  Ra tio3  as  Predicted  by  RUBY 


Explosive 

Product 

Ratio 

1.0 

Loading  Donsity,  g/cc 
1.2  1.4  1.6 

px 

TNT, 

CO i/co 

0.56 

1.00 

2.74 

7.13 

8.83 

■ti 

X 

11 

H20/Hj 

27.8 

93.9 

373 

1981 

2497 

HMX, 

COj/CO 

0.37 

0.94 

2.73 

10.52 

254 

Px  =  1.9 

Hj0/H2 

26.0 

118.4 

754 

8894 

106 

Tetryl 

COj/CO 

0.47 

1.05 

2.57 

7.76 

15.41 

px  =  1,7 

Ha0/Ha 

24.8 

91.2 

412 

2876 

104 

DATB 

COj/CO 

0.60 

1.35 

3.30 

9.38 

31.26 

px  =  1 .788 

HaO/Ha 

32.0 

119.1 

529 

3298 

104 

As  is  shewn,  the  computer  predicts  high  (Ha0/Ha)  ratios  for  C-H-N-0 

explosives  at  p  >1.0  g/cc:  so  that,  for  practical  calculational  purposes,  equi- 
0 

librium  (11)  may  be  considered  as  invariently  to  the  right  at  all  loading 
densities  under  consideration19.  The  [2  CO  C0a  +  C]  reaction,  on  the  other 
hand,  is  in  a  region  of  shifting  equilibrium  and  may  be  considered  as  predomi¬ 
nantly  to  the  right  only  at  the  higher  loading  densities  (i.  e.,  greater  than 
1.6  or  1.7  g/cc.). 

This  study  hid  as  its  original  purpose  to  develop  a  simple  method  of 
intercomparing  detonation  properties  of  experimental  C-H-N-0  high  explosives. 

Such  materials  generally  have  crystal  donsitie3  of  1.7  —  1.9  g/cc,  and  are  most 
often  used  at  high  proportions  of  theoretical  mnximum  donsity.  For  these  explo¬ 
sives  it  was  therefore  considered  a  reasonable  first  approximation  to  assume  that 
equilibrium  (10)  was  also  to  the  right  in  the  dotonation  state. 

Provided  that  no  other  species  than  Na,  H20,  CQa,  CO  and  Ha  are  present 
.in  appreciable  amounts  in  the  dotonation  gases  (and  RUDY  predicts  less  than 
2  mole  %  extraneous  Bpecies  at  1.0  g/cc,  lesser  amounts  at  higher  loading 
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cbnsiti.es),  situations  where  equilibria  (10)  and  (ll)  are  each  predominantly 
to  one  extreme  or  the  other  are  described  concisely  by  the  various  arbitrary 
assumptions  of  detonation  product  compositions.  Thus  the  [Ha0-CG-C0a]  arbitrary, 
already  mentioned,  represents  equilibrium  (10)  as  beirg  predominantly  to  the 
left  and  equilibrium  ( 1 1 )  predominantly  to  the  right;  the  [CCi-Ha0-CQa]  arbitrary 
represents  both  equilibria  as  being  predominantly  to  the  left.  In  a  like  manner, 
the  condition  in  Table  IV  at  the  higher  densities,  where  equilibria  ( 1 0)  and  (11) 
are  both  predominantly  to  the  right,  may  be  represented  by  an  [HaO-COa]  arbi¬ 
trary.  The  latter  arbitrary,  which  will  provide  the  basis  for  our  estimating 
N,  M  and  Q  in  the  present  calculations,  predicts  Na,  Ha0  and  C0a,  bit  not  CO,  as 
the  important  detonation  products,  with  Ha0  having  priority  in  formation  over  C0a, 
Given  an  explosive  compound  or  composition,  CaH.  N  0^,  in  which  there  is 
at  least  enough  oxygen  to  convert  hydrogen  to  Ha0  but  no  more  than  is  also 
req-  ired  to  convert  carbon  to  C0a,  the  [Ha0-C0a]  arbitrary  calls  for  the  formation 
of  detonation  products  according  to  the  following  decomposition  equation: 

Wc°d  — f  fN*  +  f  H.°  +  (f  -  41  co*  +  <“  '  I  +  4>  0  •  <12> 

It  follows  then  that, 


N 


arb 


2  c  +■  2  d  +  b 
48  a  +  4  b  +  56  c  +  64  d 


(13) 


M 


arb 


=  56c  ±  88  d  -  8b 
2  c  +  2  d  4-  b 


and  since. 


(14) 


-  [AH„  (detonation  products)  -  AHr  (explosive)]  \ 

o  _  _  ah  -  _ _ _ _ _ _ _ * -  , 

v*  0  formula  weight 

taking  standard  heats  of  formation  for  water  (g),  nitrogen  and  carbon  dioxide  and 
assuming  the  AH^  of  solid  carbon  to  be  nil  loads  to, 
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Q  =  «^»9b  +  47-0  (d  -  b/2)  +  AHf  (explosive) 

12  a  +  b  +  14c  +  16d 

ROBY  takes  the  NRT  term  into  account  in  its  detonation  energy  computational  for 
purposes  of  convenience,  we  ignore  it  in  the  arbitrary  calculations.  Thus,  RUBY1 8 
Q  represents  AE0|  represents  AHa.  The  difference  amounts  to  10-15  cal/g  or 
about  1%  of  Q  for  a  typical  explosive. 

Values  of  N  ^  as  calculated  from  Eq,  (13)  are  compared  in  Table  V  with 
the  corresponding  RUBY  predictions  for  a  variety  of  explosives  at  loading  densities 
from  1.0  to  1.9  g/cc. 


Table  V.  Comparison  of  Narb  with 


Explosive 

t  (rof->  "arb  . 

N 

RUBY  . 

%  dlff 

HMX 

1.903  (16)  .0338 

.0338 

0.0 

1.808  (16) 

.0338 

0.0 

1.600  (11) 

.0341 

-0.9 

1.400  (11) 

.0348 

-2.9 

1.200  (11) 

.0362 

-6.6 

1.000  (11) 

.0381 

-11.3 

PETN 

1.780  (16)  .0316 

.0318 

—0.6 

1.600  (11) 

.0320 

-1.3 

1.400  (11) 

.0333. 

-5.4 

1.200  (ll) 

.0348a 

-9.2 

Tetryl 

1.730  (16)  .0270 

.0273 

-1.1 

1.642  (16) 

.0276 

-2.2 

1.400  (11) 

.0285 

-5.5 

1.200  (11) 

.0301 

-10.3 

1.000  (11) 

.0322 

-16.1 

TNT 

1.651  (16)  .0253 

.0257 

-1.6 

1.600  (11) 

.0258 

-2.0 

1.400  (11) 

.0265 

-4.8 

1,200  (11) 

.0277 

-3.7 

1.000  (11) 

.0293 

-13.7 

DATB 

1.780  (11)  .0278 

.0279 

-0.4 

1.600  (11) 

.0281 

-1.1 

1.400  (11) 

.0287 

-3.1 

1.200  (11) 

.0298 

-6.7 

1.000  (11) 

.0313 

-11.2 

Footnote  to  Table  V: 

a)  Solid  carbon  no  longer  appears) 
the  aajne  at  lower  densities. 


NR'JBY  18  about 


NOLTR  67-66 


The  tabulated  data,  which  represent  a  fair  sampling  of  RUBY  print-out 
results,  confirm  that  N  ^  corresponds  closely  to  for  C-H-N-0  explosives 

at  higher  loading  densities,  and  differs  increasingly  from  N^ygy  as  the  densities 
decrease  and  equilibrium  (10)  shifts  to  the  left.  From  these  and  additional  data, 
average  differences  between  Narb  and  Njyjgy  ares  —0.3%  at  loading  densities  above 
1.75  g/ccj  -1.5/C  at  1.60  to  1.75  g/ccj  -3.8/6  at  1.40  g/ccj  -7.4£  at  1.20  g/ccj  and 
-11.956  at  1.00  g/cc.  Differences  between  and  Mjyjgy  and  between  and 
^RUBY  are  correspondingly  small  at  the  higher  densities  and  become  correspondingly 
greater  as  the  densities  decrease. 

Because  it  was  considered  that  the  intrinsic  inexactness  of  the  computer's 
input  information  led  to  uncertainties  of  at  least  +5%  in  RUBY's  predictions  of 
P,  we  originally  felt  that  differences  of  about  the  same  magnitude  between 
arbitrary  N,  M  and  Q  and  RUBY's  value,*?  could  be  tolerated  in  the  present  study. 

For  this  reason  the  analyses  of  Tables  IV  and  V  led  us  to  set  =1.4  g/cc  as 
a  tentative  lower  limit  of  applicability  of  the  present  calculations!,  method. 

The  results  below  shew  this  restriction  to  be  unnecessary. 

VII,  RESULTS 

Eqs.  (8)  and  (13-15)  provide  the  basis  for  a  simple  method  of  estimating 
detonation  pressures,  which  requires  is  input  information  only  the  elemental  com¬ 
position,  loading  density,  and  an  estlnata  of  the  heat  of  formation  of  the  explo¬ 
sive.  With  the  aid  of  a  desk  calculator,  a  typical  calculation  requires  less  than 
ten  minutes.  Por  comparison,  two  to  four  minutes  of  machine  time  are  required  far 
a  routine  RUBY  computation  on  the  and  a  skilled  operator  requires  an 

additional  ten  to  fifteen  minutes  to  prepare  the  input  data  and  punch  the  cards. 

Detonation  pressures  estimated  by  the  simpler  method  are  compared  in  Tahls 
VI  with  the  corresponding  RUBY  values  for  twenty  eight  materials.  Since  w«  are 
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at  this  point  concerned  only  with  reproducing  RUBY  results,  it  is  unimportant 
to  the  present  discussion  whether  input  heats  of  formation  in  the  RUBY  compute tiow> 
are  accurate;  although  the  AH^.  for  picric  acid  in  ref.  11  has  a  misplaced  decimal 
point,  we  have  used  the  same  incorrect  value  to  estimate  our  Qar-^*  Table  VI  also 
contains  several  other  instances  where  differing  estimates  of  in  refs.  11 
and  16  have  led  to  our  using  two  values  of  Q  „  for  the  same  compound. 

All  available  computer  results  for  C-H-N-0  explosives  based  on  Maderks 
parameters  and  covolumes  are  included.  In  a  number  of  cases  two  sots  of  RUBY 
results  are  listed:  those  based  on  the  RDX  parameters  and  those  obtained  using 
the  TNT  parameter  set  (Table  I).  Compounds  in  Table  VI  show  b  ranging  fran 
525  to  1728  cal/g,  M&rb  ranging  from  23.00  to  36.00  g/mole,  N&rb  ranging  from 
0.0238  to  0.0367  moles  gas/g  explosive  and  Garb  ranging  from  0.722  to  1 .000 
g  gas/g  explosive,  and  are  listed  in  order  of  decreasing  G&rb.  For  reasons  which 
will  become  obvious  from  inspection  of  the  Table  and  from  subsequent  discussions, 
the  Table  includes  results  at  loading  densities  dam  to  1.00  g/cc. 

Table  VI  near  here. 

The  hand-calculated  detonation  pressures  in  Table  VI  show  goal  agreement 
with  the  RUBY  values.  Taking  all  "uncorrected"  results,  i.  e.,  comparing  the 
hand-calculated  pressures  with  RUBY  predictions  based  on  both  parfuneter  seta,  but 
excluding  the  values  in  parentheses  for  Compounds  1-6,  differences  between 
Pcalc  8)  and  PRUBY  averaC®  +2.8956  for  the  127  data  sets  (28  explosives). 
Although  this  average  difference  falls  well  within  RUBY's  uncertainty  limits, 
the  Tablo  includes  22  date  sets  wherein  differs  from  P^y  by  more  than  % 

arxl  six  data  seta  wherein  the  difference  Is  greater  than  105(.  Since  differences 
of  comparable  magnitude f  may  also  be  observed  in  Table  V!  botween  Pjyjgy  (RBI 
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Table  VI.  Comparison  of  Pca^c  (Eq.  8)  with  for  C-H-N-0  Explosives. 

Values  of  N,  M  and  Q  Estimated  from  the  [H2Q-C02]  Arbitrary. 

■ — .  —  ■■  — . . . . .  . . - -  ■ ..  s—s—ssaau 

_  ,  ,  a  P  ,  BEK  Parameters  TNT  Parameters 

Explosive  „  calc  ,  p  ,  p 

p  (ref.)  ^arb  (Sq.  8)b  RUBY  %  diffb  RUBY  %  diff 

0  _  _  „ _ -  -  - —  -  —  -  -  -  •  -  -  ■■  — 


1.  TNM,  N  =  .0306,  M  =  32.67,  Q  =  525,  G  =  1.00C 

1.640  (11)  4.007  167.8  (157.7)  162.5  +3.3  (-3.0) 

2.  BTNEN,  N  =  .0309,  H  =32.35,  Q  =  1298,  G  =  1.000 

1.862  (16)  6.332  342.1  (321.6)  319.0  47.2  (+0.8) 

1.960  (16)  379.0  (356.3)  359.8  +5.3  (-1.0) 

3.  NG,  N  =  .0319,  M  =  31.32,  Q  =  1478,  G  =  1.000 

1.590  (11)  6.840  270.2  (254.0)  246.5  +9.6  (+3.0) 


4.  BTNEU,  N  =  .0311 ,  M  =  32.18,  Q  =  1481 ,  G  =  1 .000 

1.767  (16)  6.789  330.4  (310.6)  308.1 

1.860  (16)  365.9  (343.9  346.8 

5.  TNETB,  N  =  .0298,  M  =  33.04,  Q  =  1479,  G  =  0.985 

1.000  (16)  6.587  102,6  (96.4)  93.3 

1.200  (16)  147.8  (138.9)  133.5 

1.400  (16)  201.2  (189.1)  183.2 

1.600  (16)  262.7  (246.9)  245.5 

1.691  (16)  293.5  (275.9)  274.6 

1.780  (16)  325.2  (305.7)  306.2 

6.  PETN,  N  =  .0316,  M  =  30.41,  Q  -  1525,  0  =  0.961 

1.000  (11)  6.805  105.9  (99.5  101.6 

1.200  (11)  152.6  (143.4)  144.2 

1.400  (11  207.8  195.3)  196.4 

1.670  (11)  295.7  (278.0)  280.3 

1.770  (11  332.1  (312.2  318.8 

1.691  (16)  303.1  (284.9)  286.4 

1.780  (16)  334.3  (314.2)  321.0 

7.  KK,  N  =  .0338,  M  »  27.20,  Q  =  1481,  G  =  0.919 

1.000  (11)  6,784  105.7  107.8 

1.200  (11)  152.2  148.7 

1.400  (11)  207.2  200.5 

1.600  (11)  270.6  264.6 

1.800  (ll)  342.3  346.6 

1.712  (16)  309.9  306.0 

1.802  (16)  343.2  344.2 

8.  HMJC,  N  =  .0338,  M  =  27,20,  Q  =  1475,  G  =  0.919 

1.000  (11)  6.772  105.5  107.5 

1.200  (11)  151.9  148.3 

1.400  (11)  206.7  199.6 

1.600  (ll)  270.2  263.4 

1 .900  (11)  380.9  395.3 

1.807  (16)  344.5  346.3 

1.903  (16)  382.0  391.3 


+7.2  (+0.8) 

+5.5  (-0.9) 


+10.0  (+3.3) 
+10.7  (+4.0) 
+9.9  (+3.2) 
+7.0  (+0.6) 
+6.9  (+0.3) 
+6.2  (-0.2) 


+4.2  (—2.1) 
+5.8  (-0.6) 
+5.8  (-0.6) 
+5.5  (-0.8) 
+4.2  (-2.1) 
+5.8  (—0.6) 
+4.1  1—2.1) 


-1.9 

+2.4 

+3.3 

+2.3 

-1.2 

+1.3 

-0.3 


-1.8 

+2.4 

+3.6 

+2.5 

-3.6 

-0.5 

-2.4 


267.0 


324.0 


+10.7 


+5.6 
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Table  VI  (continued) 


9. 

DINA,  N 

=  .0333,  M  = 

27.00,  Q  =  1438,  G 

=  0,900 

1.577 

(16) 

6.561 

254.3 

247.0 

+2.9 

1.600 

(16) 

261.7 

255.1 

+2.6 

1.660 

(16) 

281.7 

276.9 

+1,7 

10. 

TNTAzB, 

N  = 

.0268,  M 

=  33.33,  <2  =  1643°, 

G  -  0.893 

1.740 

(11) 

6.271 

295.8 

300.2 

-1.5 

11.  NQ,  N  =  .0384,  M  =  23.00,  Q  =  901,  G  =  0.883 

1.691  (16)  5.528  246.3 

1.780  (16)  272.9 

12.  EDNA,  N  =  .0367,  M  =  24.00,  Q  =  1297,  G  =  0.880 

1.663  (16)  6.473  278.7  285.8 

1.750  (16)  308.9  323.0 

13.  RDX/TNT,  77/23,  N  -  .0318,  M  =  27.50,  Q  =  1436,  G  = 

1.000  (11)  6.319  93.5  100.2 

1.200  (11)  141.8  139.0 

1.400  (11)  193.0  188.1 

1.600  (11)  252.0  250.2 

1.743  (11)  299.1  304.7 

H.  HMX/TNT,  76/24,  N  =  .0318,  M  =  27.51,  Q  =  1429,  G  = 


1.000  (11)  6.305 

1.200  (11) 

1.400  (11) 

1.600  (11) 

1.809  (11) 


98.2 

141.5 

192.5 

251.5 

321.5 


99.7 

138.4 

137.4 
249.8 
333.0 


15.  DNPN,  N  =  .0322,  M  =  27.04,  W  =  1407,  G  -  0.871 

1.644  (16)  6.281  264.5  262.3 

1.730  (16)  292.9  295.1 

16.  KNsB,  H  =  .0238,  M  =  36.00,  Q  =  1728°,  G  -  0.857 

1.700  (11)  5.936  267.2  272.3 

17.  NM,  N  =  ,0319.  M  =  31.32,  Q  -  1456,  G  =  0.853 


-2.5 
— 4»4 

0.875 

-1.7 

+2.0 

+2.5 

+0.7 

-1.9 

0.875 
-1.5 
+2.2 
+2 . 6 
-KD.7 
-3.5 

+0.8 

-0.7 

-1.9 


1.128  (11)  6.769 

134.1 

130.3 

+2,9 

K/TNT,  64/36,  N  = 

.0307,  M  =  27.68, 

Q  =  1409,  G  = 

0.850 

1.000  (11)  6.063 

94.5 

95.8 

-1.5 

1.200  VII) 

136.1 

133.6 

+1.9 

1.400  (11) 

135.2 

181.3 

+2.2 

1.600  (11) 

241 .9 

240.0 

40.8 

1.715  (I1) 

277.9 

234.4 

-2.'1 

19a.  PA,  N  =  .0251,  M  =  33.06,  Q  =  1408®,  G  =  0.829 


1.000  (11)  5.415 

84.4 

87.9 

1.200  (11) 

121.4 

121.9 

1.400  (11 

165.4 

163.4 

1.600  (ll) 

216.0 

214.7 

1,760  (11) 

261.3 

265,2 

-4.0 

-0.4 

+1 .2 
+0.6 
-1.5 


248.5 

277.9 

269.4 

300.8 


-0.9 

-1.8 

+3.5 

+2.7 


288.0 


+3.7 
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Table  VI  (continued) 


19b. 

PA,  Q  = 

1261 £ 

1.672 

(16) 

5.125 

223.1 

221.2 

+0.8 

1.760 

(16) 

247.2 

243.8 

— 0.6 

20a. 

Tetryl, 

N  =  . 

0270,  M 

—  30.46,  Q  = 

I351g, 

G  =  0.822 

1.000 

(11) 

5.478 

85.3 

87.1 

-2.1 

1.200 

(11) 

122.9 

121.4 

+1 .2 

1.400 

(11) 

167.3 

1 64*6 

+1.6 

1.600 

(11) 

218.5 

218.8 

-0.1 

1.700 

(11) 

246.7 

251.5 

-1.9 

20b. 

Tetrvl, 

Q  =  1420h 

1.644 

(16) 

5.615 

236.3 

235.9 

+0.2 

1.730 

(16) 

261.8 

264.3 

-0.9 

21. 

Expl.  D, 

,  N  = 

.0285,  1 

K  =  28, 

.30,  Q  = 

:  1082, 

G  -  0.0807 

1.634 

(16) 

4.993 

207.7 

213.5 

-2.C 

207.0 

+0.3 

1.720 

06) 

230.1 

242.1 

-5.2 

229.3 

+0.3 

22a. 

DATB,  N 

=  .0278,  M  = 

28.44 

,  Q  =  12461,  G 

=  0.791 

1.000 

(11) 

5.233 

81.5 

81.7 

-0.2 

1.200 

(11) 

117.5 

114.4 

+2.7 

1.400 

(11) 

1 59.7 

155.0 

+3.0 

1.600 

(11) 

208.7 

205.8 

+1 .4 

1.788 

(11) 

260.6 

282.0 

-7.6 

264.8 

-1.6 

22b. 

DATB,  Q 

=  H51J 

1.000 

(16) 

5.030 

78.4 

76.5 

+2.4 

1.200 

(16) 

113.0 

109.9 

+2.8 

1.400 

(16) 

153.6 

151.8 

+1  »2 

1.600 

(16) 

200.3 

205.0 

-2.3 

1.745 

(16) 

238.8 

238.6 

+0.1 

1.837 

(16) 

264.3 

287.6 

—8.1 

267.8 

-1.3 

23. 

TATB,  N 

=  .0291 ,  M  = 

27.20 

,  Q  =  1075,  G 

a  0.791 

1.000 

(11) 

4.975 

77.5 

75.7 

+2.4 

1.200 

(11) 

111.7 

107.4 

+4.0 

1.400 

(11) 

152.0 

147.8 

+2.8 

1.600 

(11) 

198.5 

199.0 

-0.3 

1.895 

(11) 

278.4 

326.0 

-14.6 

297.3 

-6.4 

1.841 

(16) 

262.8 

270.1 

-2.7 

1.938 

(16) 

291.2 

304.3 

-4.3 

24. 

R-Salt, 

N  =  , 

,0345,  M 

II 

K> 

W 

• 

00,  Q  = 

1395, 

G  -  0.790 

1.520 

(16) 

6.135 

220.9 

221.2 

-0.1 

1.600 

(16) 

244.8 

246.6 

-0.7 

25. 

TNA,  N  : 

s  .0263,  M  = 

30.00, 

Q  =  1242,  G  = 

:  0.789 

1.682 

(16) 

5.077 

223.7 

219.6 

+1.9 

1.770 

(16) 

247.7 

245.6 

+0,9 

26. 

TNB,  N  : 

=  .0246,  M  = 

32.00, 

Q  =  1346,  G  = 

:  0.787 

1.604 

(16) 

5.105 

205.0 

205.4 

-0.2 

198.9 

+3.1 

1.688 

(16) 

227.0 

229.9 

-1.3 

221.1 

+2.7 
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Table  VI  (continued) 


DNPF,  N  =  .0289,  M  = 

=  26.93,  Q  =  1070,  G 

=  0.778 

1.520  (16)  4.906 

176.7 

172.5 

+2.4 

1.600  (16) 

195.8 

193.3 

+1.3 

TNT,  N  =  .0253,  M  = 

28.52,  Q  =  1282,  G  - 

:  0.722 

1.000  (k)  4.838 

75.3 

75.1 

+0.3 

76.2 

-1.2 

1.200  (k) 

108.5 

107.5 

+0.9 

106.7 

+1.7 

1.400  (k) 

147.8 

148.2 

-0.3 

145.4 

+1.6 

1.600  (k) 

193.0 

199.5 

-3.3 

194.3 

-0,7 

1.640  (11) 

202.7 

213.0 

-4.8 

205.7 

-1.5 

1.468  (16) 

162.4 

161.2 

+0.7 

1.651  (16) 

205.6 

214.4 

-4.2 

207.1 

-0.7 

Footnotes  to  Table  VI: 

a)  See  Appendix  II  for  glossary  of  compound  names  and  molecular  formulas. 

b)  Values  in  parentheses  for  Compounds  1-6  are  after  -6%  "correction"  where 
Garb  >  0,93 •  See  text. 

c)  Based  on  AHf  =  +270.4  Kcal/mole.  Probably  incorrect, 

d)  Based  cm  AHf  =  +153.8  Kcal/mole,  Probably  incorrect. 

e)  Based  on  AHp  =  -22.7  Kcal/mole.  Misplaced  decimal  in  input  data. 

f)  Based  on  AHf  =  -57.3  Kcal/mole. 

g)  Based  on  AHf  =  -15.0  Kcal/mole. 

h)  Based  on  AH^  =  +4*7  Kcal/mole. 

i)  Based  on  AHf  =  -6.0  Kcal/mole. 

j)  Based  on  AH^,  -  -29.2  Kcal/mole. 

k)  TNT  parameter  results  from  ref.  11}  RDX  parameter  results  from  ref,  16. 
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parameters)  and  PRUBy  PQrajnater8),  e.  g.,  9.7^  for  TATB  at  1.895  &/cc,  7.C# 
for  RDX  at  1.800  g/cc,  however,  it  becomes  necessary  for  each  compound  to  decide 
which  set  of  RUBY  predictions  is  more  appropriate  to  compare  with  the  results  of 
present  calculations. 

VIII.  ANALYSIS  OF  RESULTS  BASED  ON  "MORE  APPROPRIATE"  PARAMETERS 

It  has  been  mentioned  that  Mader  designed  the  STRETCH  BKW  computor  coden,ac 
so  as  best  to  reproduce  experimental  measurements  on  RDX  and  TNT.  Fundamental 
difficulties  in  finding  a  single  set  of  a,  [3,  x  and  0  in  F-q.  (4)  to  accomodate 
the  results  on  both  explosives  led  him  to  employ  two  parameters  se  ts  (Table  I), 
one  of  which  reproduced  the  RDX  measurements,  the  other  the  TNT  measurements. 

To  compute  detonation  properties  of  explosives  other  than  RDX  or  TNT,  Mader 
suggested  the  amount  of  solid  carbon  in  the  detonation  products  as  a  qualitative 
basis  for  deciding  which  parameter  set  to  use.  For  explosives  producing  greater 
amounts  of  solid  carbon  in  the  detonation,  the  TNT  parameters  were  believed  to 
be  the  more  appropriate;  for  compounds  producing  lesser  amounts  of  solid  carbon, 
the  RDX  parameters  were  considered  the  more  suitable. 

Since  the  quantity,  Garb,  in  the  present  discussions  is  an  easily  calcu¬ 
lable  measure  of  the  amount  of  solid  carbon  produced  in  the  detonation  of  an 
organic  high  explosive  (i.  e.,  weight  proportion  solid  carbon  =  1.000  -  G) ,  a 
somewhat  more  quantitative  criterion  of  parameter  suitability  suggests  itself. 

For  RDX,  G&rb  =  0.919}  for  TOT,  G&rb  =  0.722;  a  median  Gflr<b  =  0.820.  It  seems 
reasonable,  then,  that  for  explosives  with  GQrb  >  0.820  RUBY  computations  based 
on  the  RDX  parameters  might  better  accomodate  experimental  measurements,  and  that 
for  explosives  with  Garb  <  0.820  the  TOT  parameters  might  lead  to  better  corre¬ 
spondence. 

Using  Gfirb  greater  or  lesser  than  0.820  ae  a  basis  for  choosing  between 
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sot3  of  RUBY  results,  wo  may  now  compare  tho  hand-calculations  with  computer 
predictions  based  on  "more  suitable"  parameters,  i.  e.»  VG  ^RUBY 

parameters)  for  compounds  1-10  and  12-20,  vs  PRUBy  (TOT  parameters)  for 

compounds  21-28: 

Whero  Gar^  >  0.820,  19  compounds,  72  data  sets, 

Avorage  Difference  =  ^3.07%, 

Where  G  .  <  0.820,  8  compounds,  31  data  sets, 

u  Average  Difference  =  +1.77%, 

All  results,  27  compounds,  103  data  sets, 

Average  difference  =  +2.68% 

By  way  of  comparison,  differences  between  and  (l033  suitable  para¬ 

meters)  average  +3*80%  for  11  compounds,  24  data  sets. 

One  group  of  materials  deserves  special  comment.  It  may  bo  noted  that 
the  hand-calculated  detonation  pressures  are  in  all  cases  significantly  higher 
than  the  RUBY  predictions  for  the  over ox  "'zed  (Cmpds.  1-3),  C0a-balanced  (Cmpd. 
4)  and  neax^-C0a-balanced  explosives  (Cmpds.  5>6).  For  these  six  compounds  (19 
data  sets),  the  average  difference  between  Fca^c  (Eq.  8)  and  PRU0y  13  _f6»54% 
(always  positive).  In  Appendix  III  we  shall  offer  arguments  that  those  "errors" 
do  not  necessarily  reflect  any  basic  inadequacy  of  Eq.  (8) ,  but  rather  that  the 
RDX  parameters  in  RUBY  computations  may  be  inherently  unsuitable  for  noar-C0a- 
balanced  explosives. 

For  purposes  of  achieving  closer  correspondence  with  RUBY,  however  (and 
with  the  reservation  that  this  "correction"  is  not  necessarily  applicable  for  the 
prediction  of  actual  detonation  parameters),  the  following  additional  step  may 
be  incorporated  into  the  present  method  of  calculation: 

Whore  G&rb  >  0.93>  subtract  6%  from  Pcalc  (£q»  8).  (l6) 

Calculated  detonation  pressures  and  %  differences  incorporating  this  "correction" 
are  given  in  parentheses  in  Table  VI  for  Compounds  1-6.  If  these  hand-calculated 
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values  are  used,  the  results  break  down  as  follows: 

Where  G  ^  =  0.961  to  1«000,  REK  parameters  ana  a  -6%  correction, 

Average  Difference  =  +1.64% 

Where  GarB  =  Q«822  to  0.919,  EDX  parameters, 

avcrago  Difference  =  +1.82% 

Whore  G  ^  -  0.722  to  0.807,  TNT  parameters, 

Average  Difference  =  +1 .77% 

All  results,  Average  Difference  -  +1.77% 

An  alternative  method  of  comparing  the  results  of  Sq.  (8)  with  RUBY  is 
shown  in  Fig.  1,  wherein  values  of  PRUBy/?arb  arfc  Plot't;ed  against  p^.  The  solid 
line  is  of  slope  15.58  and  passes  through  the  origin,  and  thus  represents  Eq.  (8)» 
Corresponding  to  "correction"  (16),  values  of  P^ygV  are  divided  by  a  factor  of 
0.94  where  GarB  >  0.93*  It  is  seen  that  only  the  single  data  point  representing 
TATB  at  1.895  g/cc  differs  by  more  than  5%  from  the  calculated  value. 


Figure  1  rear  here 


In  the  light  of  our  ear.  ier  observation  that  Nar^,  MarB  and  correspond 
closely  to  RUBY  values  at  higher  but  not  lcwer  loading  densities  and  our  tenta¬ 
tively  having  set  p  u  1 .4  g/cc  as  a  lower  limit  of  applicability  of  the  method, 

0 

it  is  also  of  interest  to  consider  hew  differences  between  P  ,  and  PDimv  are 

calc  RUBY 

affected  try  p  .  The  breakdown  of  results,  after  the  -6%  "correction"  at  G  ,  > 

0  &rD 

0.93,  la  as  fellows: 

Average  Difference  =  +1.66%, 

Average  Difference  ■  +2.11%. 

Although  "errors"  are  indeed  slightly  greater  at  the  lower  densities,  the 
trend  is  not  strong  and  correspondence  botween  Eq.  (3)  and  RULY  is  still  sur¬ 
prisingly  good  down  to  1,00  g/cc.  That  values  of  N,  M  and  Q  from  the  [Ha0-C0a] 
arbitrary  should  differ  from  RUBY1 a  N,  H  and  Q  at  the  lower  densities  by  such 
large  amounts  as  are  shown  in  Table  V,  yet  lead  to  such  good  agreement  in  P  as 
is  shown  in  Table  VI,  raises  some  interesting  questions.  These  will  be  discussed 


Where  p  >  1 .399, 

0 

Where  p  =  1.000  to  1.200, 
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in  detail  in  the  next  paper  of  this  series.  Subsequent  papers  vdll  al^o  include 
a  comparison  of  calculations  by  the  present  method  with  experimental  detonation 
pressures  and  further  evidence  to  support  our  simplified  scheme  for  calculation 
of  detonation  velocities  (Eq.  9). 

IX.  SUMMARY  OF  THE  METHOD 

Given  the  elemental  composition,  loading  density  and  an  estimate  of  the 
heat  of  formation  of  a  C-H-N-0  explosive,  it  is  possible  to  estimate  the  detona¬ 
tion  pressure  by  the  following  simple  sequence  of  opera t;«ns: 

a)  Calculate  Narb,  Mflrb  and  fran  Eqs.  (13)»  li4)»  and  (15)»  multiply 

N  ,  by  M  ,  to  get  G  ,  . 
arb  J  arb  6  arb 

b)  Substitute  NQrb,  M&rb  and  Qarb  into  Eq.  (8)}  solve  for  P^^. 

c)  Where  correspondence  with  RUBY  (but  not  necessarily  with  actual  deto¬ 
nation  pressures)  is  desired,  subtract  6 frcn  P  if  GQrb  is  greater  than  0.93. 

At  the  current  "state-of-the-art",  such  estimates  warrant  at  least  the  same 
and  in  some  cases  possibly  greater  reliance  than  the  results  of  complex  machine 
computations  (see  Appendix  III  and  subsequent  papers). 


APPENDIX  I.  The  Dependence  of  Pressure  on  Temperature  in  the  Klstiakowsky-Wjlaon 
Equation  of  State  with  Mader's  RSK  Parameters. 

In  the  special  case  of  COa-balanced  or  overbalanced  explosives,  no  solid 
carbon  is  produced  and  G  =  1.00.  Eq.  (7),  with  Mader's  RDX  pa  race  tens,  be  canes  1 


P/N 


0.08205  T  p  [1 


+  150.1  P 

(T  +  400)0*5 


o24.0p/(T  +  40O)0,5-|  ^ 


(4-1 ) 


where  p  is  the  compressed  gas  density.  This  dependence  of  P/N  on  p  and  T  is 
shown  in  Tabls  A-I  in  the  range  of  temperatures  and  densities  of  interest  in 


detonation  calculations 
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Table  A-I.  Values  of  P/N  at  Various  Temperatures  and 
Densities 


p 

(p/nJ^sqq0 

(P/N/j-jqo0 

JO 

8 

O 

1.60 

4,099 

3,130 

2,628 

2.00 

7,025 

5,748 

5,104 

2.40 

11,270 

9,342 

9,219 

2.80 

17,245 

16,032 

15,818 

3.20 

25,468 

25,160 

26,121 

In  this  case,  N  is  relatively  independent  of  T  3irce  there  is  no  solid 

carbcn  to  participate  in  equilibria  (10)  and  (ll).  For  C~J  detonations,  P  is 

therefore  a  function  of  p  and  T  for  the  given  explosive,  the  appropriate  valuee 

being  formally  determined  by  application  of  the  hydrodynamic  equation  and  the 

C-J  condition.  However,  the  RUBY  print-outs  show  that  the  compressed  density, 

Pj,  is  mainly  a  function  of  p^  and  only  weakly  dependent  on  Q  or  the  elemental 

composition  so  that,  for  a  wide  variety  of  C-H-N-0  explosives,  RUBY's  p  's  and 

0 

pj's  satisfy  the  relationship, 

Pj/po  =  1.470  -  .05625  Pj  ,  (A-2) 

to  within  several  tenths  of  one  percent  over  the  entire  range  of  loading  densi¬ 
ties  studied. 

If  Eq.  (A-2)  is  substituted  into  (A-l)  the  result  gives  the  approximate 
P,  T  states  which  can  be  satisfied  by  a  given  loading  density.  To  show  how 
insensitive  the  detonation  pressure  ir  to  the  value  of  the  detonation  temperature 
which  is  date  mired,  extrema  values  of  T  have  been  used  to  generate  Table  A-II. 


2b 
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Table  A-II.  Dependence  of  P  on  T  in  Eq.  {A— 1 )  at  Various 
Loading  Densities. 


p 

0 

P4500O//p2100° 

P4500o//pl200° 

1  .’6 

1.310 

1.560 

1.47 

1 .222 

1.376 

1.80 

1.145 

1 .222 

2.13 

1.076 

1.090 

2.«,8 

1.012 

0.975 

Since  most  detonation  temperatures  calculated  by  RUBY  fail,  within  the 
range  2000  to  4000°K,  it  is  apparent  frcm  Table  A- 1  that  an  estimate  of  the 
detonation  temperature  to  within  10^  at  the  lower  densities  would  fix  Pj  to 
within  about  ±2-3%  and  at  higher  densities  the  error  would  became  progressively 
smaller* 

The  more  general  case  where  G  /  1.00  requires  a  more  complex  analysis  of 

interrelationships  between  RUBY’s  p  ,  pT,  p  ,  N  and  V  and  between  P,  N,  G  and 

0  w  g  s  s 

T,  but  the  conclusions  are  essentially  the  same.  Dependence  of  P  ou  T  is  not 

strong  at  the  lower  densities  and  becomes  progressively  weaker  as  p  increases, 

0 

until  at  p  =  ca  2.4  g/cc  the  inversion  in  sign  of  AP/AT  is  observed.  Such  very 
0 

weak  dependence  of  P  on  T  above  1.80  g/cc  as  is  shown  in  Table  A-II  may  signal 
a  lessening  adequacy  of  Mader’s  RHX  parameters  at  the  higher  densities  and  account 
for  the  upward  trend  in  Ppygy  relative  to  (Eq.  8)  at  1.85-1.96  g/cc  in 

Figure  1 . 

Such  results  may  also  suggest  that,  although  excellent  as  an  interpolation 
tool  as  was  intended  by  RUBY's  designers  (i.  a.,  for  explosives  with  compositions 
and  properties  between  REIX  and  TNT,  see  also  Appondix  III),  the  RUBY  code  may  be 
less  satisfactory  in  extrapolative  situations.  Various  RUBY  users  have  computed 

P7 
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detonation  properties  of  hypothetical  explosives  at  predicted  densities  as  high 

as  2*1-2. 2  g/cc,  and  have  used  the  resell ts  as  a  basis  for  extended  synthesis 

programs.  It  is  now  suggested  that  predictions  of  explosive  properties  based 

on  such  computations  are  subject  ti  serious  question* 

Similar  analyses  of  interrelationships  between  other  of  the  quantities  in 

RUBY  print- outs  which  interact  to  produce  Pj^y  have  been  carried  out.,  These 

have  allowed  a  series  of  approximations  whereby  Eqa.  (7)  and  ( A— 1 )  are  modified 

to  yield  still  another  expression,  which  has  its  roots  in  the  behavior  of  the 

K-W  equation  of  state,  but  which  closely  reproduces  the  P-N-M-Q~p  relationships 

0 

in  empirical  Eq.  (8),  The  reasonorg  behind  these  approximations  is  rather 
involved  and  c?  probable  interest  to  only  a  limited  group  of  readers}  thair 
detailed  discussion  will  therefore  be  deferred  to  a  subsequent  paper  in  this 
series. 

APPENDIX  II.  Glossary  of  Compound  Names  and  Molecular  Formulas* 

1.  TNM,  CN40a,  tetranitr one thane 

2.  BTNEN,  bis(2,2,2-udrdtroethyl)nitraraine 

3.  NG,  C3H3N309,  nitroglycerine,  glycorol  trinitrate 
4*  BThEU,  C5^Na013,  bi3(?,2,2-trlnitroethyl)urea 

5*  TVETB,  C6HaN6014,  2,2,2-trinitroetliyl  4,4,4-trinitrobutymte 

6.  PETN,  C3HaIJACia,  pentaorithritol  tetranitrate 

7.  RDX,  C3H6N60iS,  cyclotriniothylone  trinitramino,  1 ,3,5-triaza-1 ,3,5-trinitro- 

cyclohexane 

8.  HMX,  C4HaNa0a,  cyciototrQ.iiothvlQno  te  tra  nit  rami  ne,  1 ,3,5»7-tetrrt»a-1 ,3»5*7- 

totranitrocyclooc  tane 

9.  DINA,  CAHaNAOa,  di(2-nltroxyet-hyl)nitramine 

10.  TNTAzB,  C6Nia06,  1 ,3,5-tria-/,ido-2,4»6>-trinitrobenzeno 

11.  NQ,  CH^N^Oa,  nitrogvanidine 


NOI.TR  67-66 


12.  EDNA,  CaH6N404,  a thy lone  dinitramine,  1 ,2-di( nitramino) ethane 

11.  RDX/TNT,  77/23,  C3 .04H7 .46N6. afl07 ,73,  Cylcotol 

14*  HMX/TNT,  76/24,  0.03^9. 0.43 »  Gctol 

15*  DNPN,  C^H^qN^O^q,  bia{2,2— dinitropropyl)nitramine 

l6.  HNsB,  CjH^O^,  haxanitrosobanzene,  benzotrifuroxan 

17'.  NM,  CHjNOj,  nitrametbane 

1H.  RDX/TOT,  64//36,  C6,a3Ha#73N7j6309,30#  Composition  B 

19.  PA,  C^HjNjO.,,  picric  acid,  2,4,6-trinitrophenol 

<:0.  Tetryl,  C7H3N308,  N-methyl-N-nitro-2,4,6-trinitroaniLina 

21.  Expl.  D ,,  C6K*N*07 ,  ammonium  picmte 

.^2.  DATB,  CjH3N306,  1 ,3-diamino-2, 4, <>■  trinitrobenzene 

'•3»  TATB,  C6H^N606,  1 ,3,5-t:rlaminv>-2,4,6“trinitrobenzene 

24.  It-Salt,  C3 HgN60j ,  cyclotrimethylene  trinitroaamine,  1 ,3,5-triaza-1 ,3,5- 

trinitrosocyclohexane 

25.  TNA,  C6H4Na06,  2,4,6"trinitroanilin0,  picrajnide 

26.  TN3,  C6H3N306,  1 ,3,5-trinitrobenzane 

27.  DNPF,  CjoHja^Oia*  bis(2,2-d.' aitropr'opyl)  fumerate 

28.  TNT,  07^11,0$,  2, 4, 6- trinitrotoluene 


APPENDIX  III.  Tho  -6%  "Correction"  at  Garfr  >  0.93. 

The  reasoning  behind  our  suggestion  that  the  need  for  the  -6£  correction 
for  heax^COj- balanced  explosives  is  not  necossarily  in  consequence  of  any  basic 
inadequacy  of  Eq«  (8)  is  as  follows:  At  most  loading  densities  currently  under 
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consideration,  RGX  parameter?  lead  to  higher  values  of  than  TNT  parameters 

(see  Table  VI),  the  difference  increasing  at  the  higher  loading  densities21. 

It  has  been  mentioned  that,  extending  Mader'e  reasoning,  REK  parameters  are 
"ax.ctly  suitable"  at  G  ^  =  0,919  and  TNT  parameters  are  "exactly  suitable"  at 
Garb  “  °  .722.  At  other  values  of  G  ^  neither  parameter  set  is  "exactly  BuitabLe" 
but,  setting  Garb  greater  or  lesser  than  0.820  as  a  criterion  for  choice,  one  or 
the  other  parameter  set  is  "more  appropriate". 

"Exact  suitability"  at  all  values  of  G  ^  would  require  either  an  infinite 
number  of  K-W  parameter  sots,  or  equations  which  adequately  expressed  the  para¬ 
meters  as  continuous  functions  of  G  or  some  other  appropriate  property  of  the 

explosive.  Although  we  are  not  now  in  a  position  to  offer  such  equations,  it 

# 

is  nevertheless  profitable  to  consider  a  quantity,  which  would  represent 

the  pressure  predicted  by  the  computor  if  given  as  input  information  "exactly 
suitable"  values  of  a,  j3,  v  and  9. 

From  the  relationship  between  RDX -parameter  RUBY  results  and  TNT-p»rameter 


resultsj  it  follows  that  at  increasing  values  of  Gar^  between  0.722  and  0.919» 
^RUBX  parameters)  should  3how  a  decreasingly  positive  bia3,  and  (TNT 

parameter)  should  show  an  increasingly  negative  bias  relative  u>  Ppyjgy8  On 
detailed  examination  of  Tabic  VI,  such  trends  between  Ppggy  and  P  ^  (Eq.  8)  may 


readily  be  discerned  and  it  is  now  suggested  that  average  differences  would  be 

« 

even  smaller  if  carnpjarisons  wore  between  P  &nd  such  a  Ppjgy, 


Extending  the  same  reasoning,  P^^y  (HDX  parameters)  should  run  increasingly 
* 

low  relative  to  at  increasing  values  of  Garb  abovo  0.919.  In  other  words, 


givon 


sultablo"  parameters,  the  computor  would  predict 


pressures 


for  overotxldized,  C02-- balanced  and  neai-CO^-balanced  explosives  than  are  reported 
in  refs,  11  and  16.  Although  we  cannot  now  say  whether  such  increasos  would  be 
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sufficient  to  offset  the  "errors"  in  Table  VT  for  compounds  1-6,  we  strongly 
auspect.  that  such  might  be  the  case. 
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Footnotes  and  References  (continued): 

19)  The  position  of  this  equilibrium  being  a  sensitive  function  of  density,  how¬ 
ever,  it  would  tend  to  introduce  complications  at  p  <1.0  g/cc. 

20)  C,  L.  Mader,  Detonation  Performance  Calcula tlons  UsinR  the  Kistlakowsky- 
Wilson  Equati on  of  State.  Los  Alamos  Scientific  Laboratory  Report,  LA-2613 
(1961)'. 

21 )  Differences  between  RUBY  computations  based  on  the  two  parameter  sets  depend 
strongly  on  leading  density.  At  about  1.15  g/cc,  both  parameter  sets  give 
about  the  same  results.  As  loading  densities  increase  above  this  value,  the 
RDX  parameters  give  increasingly  higher  values  of  P  than  the  TNT  parameters# 
and  at  lower  densities  the  converse  is  the  case  but  to  a  lesser  extent. 
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